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CLAIMS 



I claim: 

A method for seismic data processing, comprising: 

constructing explicit depth extrapolation operators with variable operator lengths 
depending on maximum dip angle, accuracy condition, and wavenumber; and 
constructing operator tables using the explicit depth extrapolation operators. 

The method of claim 1, further comprising: 

performing depth migration using the explicit depth extrapolation operators from the 
operator tables. 

The method of claim 1, wherein the step of constructing operator tables comprises: 
selecting a maximum dip angle; 
selecting a type of operator; 

selecting an accuracy condition for the type of operator; 
selecting a plurality of wavenumbers; and 

performing the following steps for each of the plurality of wavenumbers: 
selecting a plurality of operator lengths for the selected wavenumber; and 
performing the following steps for each of the plurality of operator lengths: 

determining if the selected type of operator with the selected operator length 
satisfies the selected accuracy condition at the selected maximum dip angle at the 
selected wavenumber; 

determining if the operator length is the smallest operator length satisfying the 
accuracy condition at the selected maximum dip angle for the selected 
wavenumber; and 

storing the operator length in an operator table if the operator length is the 
smallest operator length satisfying the accuracy condition at the maximum dip 
angle for the selected wavenumber. 
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The method of claim 3, wherein the steps performed for each wavenumber comprise: 
selecting a plurality of pairs of first operator half-lengths and second operator half- 
lengths for the selected wavenumber; and 

performing the following steps for each of the plurality of pairs of first operator half- 

lengths and second operator half-lengths: 

determining if the pair of first operator half-length and second operator half-length 
satisfies the selected accuracy condition at the selected maximum dip angle at the 
selected wavenumber; 

determining if the selected type of operator with the first operator half-length is the 
smallest first operator half-length satisfying the accuracy condition for the selected 
wavenumber and the second operator half-length is the smallest second operator half- 
length satisfying the accuracy condition for the selected wavenumber; and 
storing the pair of first operator half-length and second operator half-length in an 
operator table if the first operator half-length is the smallest first operator half-length 
satisfying the accuracy condition for the selected wavenumber and the second 
operator half-length is the smallest second operator half-length satisfying the 
accuracy condition at the selected maximum dip angle for the selected wavenumber. 

The method of claim 1, wherein the step of constructing operator tables comprises: 
selecting a maximum dip angle; 
selecting a type of operator; 

selecting an accuracy condition for the type of operator; 
determining preliminary operator tables for varying operator lengths; 
selecting a plurality of wavenumbers; and 

performing the following steps for each of the plurality of wavenumbers: 

selecting a smallest operator length in the preliminary operator tables satisfying the 
accuracy condition for the selected wavenumber; and 
storing the smallest operator length in an operator table. 

The method of claim 4, wherein the steps performed for each wavenumber comprise: 
selecting a pair of smallest first operator half-length and smallest second operator half- 
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7. 



length in the preliminary operator tables satisfying the accuracy condition for the selected 
wavenumber; and 

storing the pair of smallest first operator half-length and smallest second operator half- 
length in an operator table. 

The method of claim 2, wherein the step of performing depth migration comprises: 
selecting a seismic data set; 

selecting a plurality of depths in the seismic data set; 
selecting sub-domains for each of the plurality of depths; 
determining a velocity model for the seismic data set and the plurality of depths; 
determining a lowest velocity for each sub-domain from the velocity model; 
selecting an operator table with varying operator lengths; 
selecting a plurality of frequencies in the selected sub-domain; and 
performing the following steps for each of the plurality of frequencies: 
performing the following steps for each of the plurality of depths: 

performing the following steps for each of the sub-domains at the depth: 

calculating a highest wavenumber in the sub-domain from the frequency and 
the lowest velocity for the sub-domain; 

selecting a maximum operator length for the calculated highest wavenumber 
from the operator table; 

applying operators from the operator table with operator lengths no greater 
than the maximum operator length to the selected sub-domain. 

The method of claim 7, wherein the operator table with varying operator lengths 
comprises explicit extrapolation operators with varying operator lengths. 



8. 



9. The method of claim 7, wherein the step of selecting an operator table with varying 
operator lengths further comprises: 

interpolating the operator table with varying operator lengths. 



10. 



The method of claim 7, wherein each operator length comprises a pair of first operator 
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half-length and second operator half-length. 

11. The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 

P(x i ,y j ,o),z + Az) = 

where P(x„ # ft), z+Azj is a seismic wavefield at lateral location x, = iAx, # =/Av, and depth 
z+Az; co is angular frequency; P(x h y h co, z) is a seismic wavefield in a sub-domain D at depth z; 
Ax, Av, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
L(k(o max (D, z)) and M(k w max (D, z)) are operator half-lengths in the x- and y- coordinate directions, 
respectively, for the highest wavenumber k m max (D, z) in the sub-domain D at depth z; W(x u yj, 
kjx,, yj, z), Az) is an explicit depth extrapolation operator; and hjxu yj, z) is a local wavenumber. 

12. The method of claim 1 1, wherein the wavenumber kjxi, y jt z) is given by the formula: 

. ft) 

*•<*<• ^ 2) = 7^^)' 

where c(x it y jt z) is a local propagation velocity. 

13. The method of claim 11, wherein the highest wavenumber KT(D, z) for a frequency co 
for a sub-domain D at depth z is given by the formula: 

where c min (D, z) is the lowest velocity for the sub-domain D at depth z. 



26 



PGS-03-03US 



14. The method of claim 3, wherein the step of step of selecting a type of operator comprises 
applying the following formula: 

/>(*,., <y,z + Az) = X W(x l ,k e) (x i ,z),Az)P(x i -x l ,co,z), 
;=-i(C x (Cz)) 

where P(x b co, z+Az) is a seismic wavefield at lateral location x, = iAx and depth z+Az, co is 
angular frequency, P(x h co, z) is a seismic wavefield in a sub-domain D at depth z, Ax and Az are 
step lengths in the x- and z-coordinate directions, respectively, L(k w max (D, z)) is an operator half- 
length in the x-coordinate direction for the highest wavenumber k a max (D, z) in the sub-domain D 
at depth z; W(x it h,(xi, z), Az) is an explicit depth extrapolation operator, and k a ( Xi , z) is a local 
wavenumber. 

15. The method of claim 14, wherein the wavenumber k m (x u z) is given by the formula: 

co 

where c(x it z) is a local propagation velocity. 

16. The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 



/>(*,. ,y.,a>,z + Az)= X W > (*• ( x < ' yj ' z)> ^ h < (x ' ' y J t(0 ' z) ' 

/=0 



where P(x u y h co, z+Az) is a seismic wavefield at lateral location x, = i Ax, yj =j Ay, and depth 
z+Az; co is angular frequency; P(x u y h co, z) is a seismic wavefield in a sub-domain D at depth z; 
Ax, Ay, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
Kk^fD, z)) is an operator length for the highest wavenumber K max (D, z) in the sub-domain D 
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at depth z; W\(ka>(xu yj, z), Az) is an explicit depth extrapolation operator expressed in cylindrical 
coordinates; k m (x it y Jf z) is a local wavenumber; and hi(x if yj> (o y z) are auxiliary fields resulting 
from applying recursive 2D Chebyshev filters to the seismic wavefield P(x it y Jt co, z). 
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